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Jet finding is a type of optimization problem, where hadrons from a high-energy collision event 
are grouped into jets based on a clustering criterion. As three interesting examples, one can form 
a jet cluster that (1) optimizes the overall jet four-vector, (2) optimizes the jet axis, or (3) aligns 
the jet axis with the jet four-vector. In this paper, we show that these three approaches to jet 
finding, despite being philosophically quite different, can be regarded as descendants of a mother 
optimization problem. For the special case of finding a single cone jet of fixed opening angle, the 
three approaches are genuinely identical when defined appropriately, and the result is a stable cone 
jet with the largest value of a quantity J. This relationship is only approximate for cone jets in the 
rapidity-azimuth plane, as used at the Large Hadron Collider, though the differences are mild for 
small radius jets. 


I. INTRODUCTION 

Jet algorithms are essential tools for connecting long¬ 
distance measurements made on hadrons to short- 
distance interpretations based on perturbative quarks 
and gluons. Since there is a fundamental mismatch be¬ 
tween color-singlet hadrons and color-carrying partons, 
there is no way to define an ideal jet finding procedure. 
For this reason, a variety of jet algorithms have been in¬ 
troduced with different underMng philosophies and dif¬ 
ferent practical advantages [U, 0 . 

In this paper, we expose a surprising connection be¬ 
tween three seemingly unrelated ^proaches to jet find¬ 
ing: jet function maximization l-iettiness mini¬ 
mization and stable cone finding |10l - [T^ . all re¬ 

viewed in Sec. [TTl Philosophically, these algorithms are 
quite different, so one might think that they would yield 
rather different jets. Instead, all three algorithms yield 
(approximately) conical jets of (approximately) fixed ra¬ 
dius R, with a high degree of correlation between the 
methods. As we will show, this correlation is not an ac¬ 
cident, since jet function maximization and 1-jettiness 
minimization can be viewed as descendants of a mother 
optimization problem, whose solution is a stable cone jet. 

This relationship is most transparent in electron- 
positron collisions, where jets are typically defined in 
terms of particle energies and angles. Remarkably, with 
appropriate definitions, all three algorithms can be made 
to yield identical cone jets. We will prove this in Sec. Hill 
by showing how these three methods can be derived from 
optimizing a common meta function. This optimization 
is equivalent to finding all stable cones and then choosing 
the one with the largest value of J (defined below). In 
Sec. lYl we explain this relationship in more detail by 
performing a two particle case study. 

Turning to the Large Hadron Collider (LHC), the in¬ 
trinsic differences between the three algorithms become 
apparent. In proton-proton collisions, jets are typically 
defined in terms of particle transverse momenta and 
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rapidity-azimuth distances. As discussed in Sec. |Vl the 
transverse momentum of a jet is not exactly the same as 
the summed transverse momenta of its constituents. For 
this reason, there is now a difference between optimizing 
a four-vector (as in jet functions), optimizing a light-like 
axis (as in 1-jettiness), and aligning the jet axis with the 
jet momentum (as in stable cones). Despite these dif¬ 
ferences, though, we show that the algorithms still give 
similar results when the jet radius R is small. 

II. REVIEW OF JET ALGORITHMS 

To begin, we briefly review these three approaches to 
jet finding, all of which are infrared and collinear (IRC) 
safe. For simplicity, in this paper we only consider finding 
the single hardest jet in an event, though all of these 
methods can be adapted to identify multiple jets. For 
consistency of notation, we always use R to refer to the 
adjustable jet radius parameter in each algorithm, which 
can be identified with the true jet radius only in the small 
R limit. 

• Jet function maximization il. Here, the jet find¬ 
ing strategy is to find the subset of particles in an 
event that maximizes a jet function J(P^), where 
Pf^ is the four-vector of the candidate jetQ The 
original paper Q introduced a jet function appro¬ 
priate for electron-positron collisions: 

1 777 2 

Joe.,{Pu) = E--^—, ( 1 ) 

where E and m are the total energy and mass of the 
candidate jet, and the original paper used the nota¬ 
tion R^ -T 1//3. Maximizing this Jorig gives quasi- 
conical jets with nearly fixed radius ~ R. Since 
then, this algorithm has been adapted to collisions 


^ The name “jet function” and the symbol J should not be con¬ 
fused with earlier usage in the context of factorization and re¬ 
summation, e.g. lUllSlI . 
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at the LHC @ to yield nearly conical jets in the 
rapidity-azimuth plane. 

• 1-jettiness minimization This approach is 

based on finding a light-like jet axis n = (1, h) that 
minimizes 1-jettiness 7i(n). For electron-positron 
collisions, one possible choice of 1-jettiness measure 
is (see also i) 

Ti{n)= min|F;,,^^^| , (2) 

where pi = (Ei,pi) is the four-vector of particle 
i, and one often sees the notation — ?> 2p. The 

minimum inside of 7i partitions the event into an 
unclustered (or “beam”) region and a jet region, 
and after minimizing over n, the jet region is ex¬ 
actly conical (for massless particles) with radius 
~ R. This approach has been generalized to N- 
jettiness jet finding through the XCONE jet algo¬ 
rithm [3,[I3- 

• Stable cone finding [Tol - [l^ . Algorithms like SIS- 
CONE [T^ search for all conical jet regions of radius 
R which are stable, meaning that the total jet mo¬ 
mentum is aligned with the jet axis. This algorithm 
works equally well for opening angle as for rapidity- 
azimuth distance. From the set of stable cones, 
one must then choose the desired jet(s). By con¬ 
struction, the hardest jet after pro gres sive removal 
(e.g. SISCone-PR as described in [1^) is a perfect 
cone. While “hardest” typically refers to the jet 
with highest (transverse) momentum, we will find 
it interest]^ to consider the jet with largest J. As 
shown in [ll|| , stable cone finding can be viewed as 
an optimization problem, which is closely related 
to 1-jettiness minimization Q. While historically 
there were issues with IRC safety of SMcific itera¬ 
tive cone algorithms (see discussion in [l|,0), those 
issues can be resolved using seedless methods (like 
in SISCone) or IRC safe seeds (like in XCone). 


III. A META OPTIMIZATION PROBLEM 


where R will soon be identified with the same parameter 
in ©■ Here, we work with massless final-state parti¬ 
cles with Ei = \pi\] one can adapt this construction to 
massive particles by performing measurements in, say, 
the p-scheme or A-scheme where massive particles are 
replaced by massless proxies (see e.g. [la,[^). We claim 
that maximizing M over all possible and n simul¬ 
taneously yields a jet function maximum, a 1-jettiness 
minimum, and a stable cone jet. 

First, consider maximizing M with respect to n, keep¬ 
ing the candidate jet fixed. Introducing a Lagrange 
multiplier A to enforce unit norm 

M(P^,n) ^M(P^, n) + A(n2-l), (4) 

one can show that M is maximized by 

n°P‘=(l,-^). (5) 

\ \P\) 

Note that n lives on a smooth compact space (i.e. the 
surface of a unit sphere), so there are no edge conditions 
to check. Plugging this back into ([3]), we find 

M(P^,n°P‘) = P - = J(P^). (6) 

Thus, maximizing M over {P^,n} is equivalent to maxi¬ 
mizing the jet function J over P^. Since 

2(P-|P|):^!^ (7) 

for m E, this jet function is very similar to the original 
definition in ([T|) □ Following the logic of Q , a non-trivial 
J maximum is guaranteed to exist since there are a finite 
number of P^ partitions and a single particle has a larger 
J value than the empty partition. 

Next, consider maximizing M over the candidate jet 
P^i, keeping n fixed. Crucially, M is linear in the final 
state particles, so we can write M as 

M(P^,n)= (9) 

zGjet 


Since the three above approaches yield (nearly) coni¬ 
cal jets and all are based on optimization, it is perhaps 
not surprising that, with suitable modifications, they can 
yield identical jets. Here, we show this explicitly for 
cone jets of fixed opening angle, as relevant for electron- 
positron collisions. We discuss obstructions to generaliz¬ 
ing this to the LHC in the Sec. |Vl 

Consider the following meta function which depends 
both on a candidate jet P^ = (P, P) as well as on an 
auxiliary light-like axis n = (l,h): 


where “jet” refers to the set of particles contributing to 
Pn (i-e- EiejetP* = P)- ^ny final-state 
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Ref. introduced a generalization of Q with an extra param¬ 
eter n, 


j(") = £" ( 1 - 


n m? A 

£2 j 


(8) 


where jn —>■ 1//3 in the original notation. Because jd is 
not linear in the jet kinematics, we know of no meta function 
that can directly handle this case. That said, maximizing is 
equivalent to maximizing V , and J in the m E 

limit. 


M{P^,n)=E 


P2 ’ 
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particle that contributes positively to M will increase 
M’s value, so for the optimal jet we must have 

M{P°^\n)= ^ maxji;, - , (10) 

2 G event ^ 


where now the sum runs over all particles in the event. 
Rearranging this formula as 

E - M{P°P\n) = min|£;„^^^| =ri(n), 

( 11 ) 

we recover precisely ©. Thus, taking into account the 
minus sign, maximizing M over {P^,n} is equivalent to 
minimizing 1 -jettiness 71 over n. 

Finally, we can maximize M with respect to both n 
and by combining these two analyses. Since M has 
a maximum, there must be an optimal subset (per¬ 
haps more than one in degenerate phase space configu¬ 
rations) and a corresponding optimal axis n°P*. Fixing 
P°E and minimizing 7 i(n) with respect to n using the 
same Lagrange multiplier trick in Q, the optimal axis 
must be 


i°E = 1 


popt 

|Popt| 


( 12 ) 


Now looking at (TTOl) with fixed n°P*, the optimal subset 
P°P* must consist of all particles with 


2 n°Pt • p, 
E, 


< p2. 


(13) 


This is precisely the condition for a stable cone. To see 
why it is stable, note that (TT^ implies that the jet axis 
n°P* is aligned with the jet momentum P°p* . To see why 
it is a cone, note that for massless particles with \pi\ = Ei, 


2n-pi 

E, 


2 2 cos j 


(14) 


where 9n,i is the angle to the axis. Thus, the condition in 
(IT3l) specifies particles contained in a cone of half opening 
angle Ptrue around n°P* with the identification 


\/2 - 2 cos Ptrue = P, (15) 


where Ptrue — P at small radii. 

We have therefore demonstrated that optimizing 
M(P^,n) is equivalent to optimizing the following two 
functions: 


j(Pm) = e- 


2{E-\P\) 


P2 


Ti(n) = Y min<^Pz, 


zGevent 


2n-pi 

P2 


(16) 

(17) 


Moreover, the optimal jet P))p* is specified by a stable 
cone centered on the axis h°P'’ oc P°p‘. Of course, ac¬ 
tually finding the optimal conhguration is a challenging 


computational problem. In practice, one can make use 
of the converse statement that among all stable cones, 
the one that has the largest value of J(P^) is guaranteed 
to be the one that optimizes MfP,,, n). Therefore, one 
can run spherical SISCone-PR [i2,[3 1° stable 

cone jets, and then simply select the jet with the largest 
value of J{P^i)- 


IV. TWO PARTICLE CASE STUDY 

To better understand the above analysis, it is instruc¬ 
tive to study the simplest case of two massless particles 
with equal energies P/2, separated by opening angle 6 * 12 . 
For ease of notation, we define 

y 2 - 2 cos^ = v' 2 - 2 cos 0 i 2 = Pi 2 , (18) 

such that 012 > P 12 > P 12 , but all three agree in the 
small angle limit. 

In the jet function approach using (HU), there are three 
candidate jets to consider: particle 1 alone, particle 2 
alone, or particles 1 and 2 together. The corresponding 
jet function values are 

JiK) = JiPi'i = I. (19) 

+ = ( 20 ) 

Choosing the maximum J value, the condition for the 
two particles to be clustered into a single jet is 

J maximum: R 12 < V2R. (21) 

In the 1-jettiness approach, it is straightforward to 
prove that the axis that minimizes 7i(n) must lie along 
the line between particles 1 and 2 , such that 


0rt,l + 0rt,2 — 012- (22) 

Furthermore, one can show that local minima of 7i(n) 
can only appear at three possible points: 

022 , 1 = 0 , 0 ft. 2 =O, 9^,1 =0n,2 = ^, (23) 

corresponding to the three candidate jet configurations 
above. For R 12 < R (note the tilde), only 9n,i = 9^:2 is a 
true local minimum, so the two particles are always clus¬ 
tered together. For R 12 > 2R (no tilde), only 9n^\ = 0 
and 0 ft ^2 = 0 are true local minima, so the two parti¬ 
cles are never clustered together. For intermediate i?i 2 
values, all three minima are present with corresponding 
1 -jettiness values of 

Ti (0ft,i=0) = ri (022.2 =0) = |, (24) 

ri(0ft.i=0ft.2) = |^. (25) 
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Choosing the global 7i minimum, the particles will be 
clustered if 

global 7i minimum: R 12 < V2R, (26) 

yielding the same result as the jet function approach. Al¬ 
ternatively, one might be content with only finding a local 
7i minimum. For example, XCone [l^, [l^ searches for 
local 7i minima using IRC safe seeds based on fc^-style 
clustering [U, , leading to the clustering condition 

fer-seeded local 7i minimum: R 12 < R. (27) 

With IRC safe seeds, this is a perfectly acceptable jet 
finding strategy from the point of view of perturbative 
QCD calculations 0 

Turning to the stable cone approach, the three local 
minima in 7i correspond precisely to the three possi¬ 
ble stable cone configurations. The default behavior of 
spherical SISCone-PR [H, [T^ is to progressively re¬ 
move the stable cone with the largest energy, meaning 
that the two particles will always be clustered if they fall 
within a common cone: 

A-ordered SISCone-PR: i?i 2 < 2R. (28) 

To mimic precisely the behavior of J maximization or 7) 
minimization, though, one has to progressively remove 
the jet with the largest value of 7(P^): 

J-ordered SISCone-PR: R 12 < V2R. (29) 

Only for this version of stable cone finding does one re¬ 
cover the remarkable equivalence derived in Sec. imi 

As a historical note, the presence of local 7i minima 
led to many Tevatron-era discussions about cone jet al¬ 
gorithms [lol . [H, where a variant of 7i was known as 
the “Snowmass potential”. One issue is that of midpoint 
seeds, since if one only uses the location of the two parti¬ 
cles themselves as seeds, then one misses out on a possible 
stable cone centered on their midpoint, as in (I23p . Since 
the location of 71 minima is an IRC safe property, this is 
the reason why seedless methods are preferred. A related 
midpoint issue is that of disappearing minima, since in 
data, where the two particle directions are smeared out 
by showering and hadronization effects, the local 7i mini¬ 
mum at the midpoint disappears for R 12 ^ RsepR, where 
Rsep — 1.3 [ 2 ^. Disappearing minima do not present an 
issue for IRC safety, but they do lead to poor conver¬ 
gence of perturbation theory. In that context, it is in¬ 
teresting to note that 1.3 is comparable to the factor 
in (1^ . This suggests that 7-ordered cone finding may 
be more robust to showering and hadronization effects 
than A-ordered cone finding, since 7-ordering aims to 
find global 7i minimum instead of just a local one. We 
leave a detailed analysis of this possibility to future work. 


® Crucial to achieving IRC safety for XCONE, the number of seeds 
is always equal to the desired number of jets 'N, independent of 
the final state particle multiplicity. 


V. DIFFERENCES AT THE LHC 


For LHC applications, it is typically advantageous to 
work with transverse momenta and rapidity-azimuth dis¬ 
tances. Here, jet function maximization, 1-jettiness mini¬ 
mization, and stable cone finding are indeed inequivalent, 
which we now show. 

Jet functions are based on optimizing the jet four- 
vector as a whole. Therefore, we expect the jet function 
to be a function of the total transverse momentum Hp, 
total mass m, or the combination Et = \/P^^+rn?. For 
example, the chosen jet function in [g was 

1 TTl^ 

A.d{P^) = ET-j^ — , (30) 

where i? —> 1/in the original notation. As explained 
in 0, while the jet region is nearly conical when using 
7had, the jet momentum direction is somewhat offset from 
the jet center at forward rapidities. 

By contrast, 1-jettiness is based on optimizing a light¬ 
like axis, using a measure that is additive over the indi¬ 
vidual particles i. As shown in Q (see also dill), a 
cone algorithm can be defined via 


pnininjl,^^! , (31) 

2 G event ^ '' 

where R is the jet radius and j3 is an angular exponent 
(not to be confused with /3 in the jet function literature). 
For any value of /3, the jet region is perfectly conical and 
the jet axis is located at the center of the jet region. For 
the special value of /3 = 2, the jet axis is aligned with the 
PT"weighted jet momentum (see (1551) below). 

At first glance, one might think that Jhad and 
could be massaged into a common meta function Mhad 
by just making suitable E ^ pT replacements. However, 
there is a fundamental mismatch due to the fact that 


For this reason one cannot construct a suitable iWhad 
that depends solely on and n, since one also needs 
information about the scalar px sum within the jet. Of 
course, one could find an Mhad if one relaxes the require¬ 
ments that 7had depends solely on or that T^^ad 
be expressed as a sum over i, though that goes somewhat 
against the original philosophies of those approaches. 

For stable cone finding, there are two different defini¬ 
tions in common use. In m , a stable cone is defined 
such that the jet axis aligns with the i?T-weighted cen¬ 
troid of its constituents. For massless constituents, this 
is the same as minimizing in SB) with /3 = 2 [^, 
yielding a jet axis located at rapidity/azimuth location 


y = 


SiGjet Vi PTi 


Sigjet PTi 
^zCjet PTi 


(33) 
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J E'J' 

XCONE N = 1 

SISCone-PR 

anti-kt 

C/A 

kt 

J E'j' 

— 

93% 

96% 

93% 

81% 

74% 

XCONE N = 1 

93% 

— 

94% 

96% 

85% 

77% 

SISCone-PR 

96% 

94% 

— 

94% 

82% 

75% 

anti-fct 

93% 

96% 

94% 

— 

84% 

76% 

C/A 

81% 

85% 

82% 

84% 

— 

83% 

kt 

75% 

78% 

76% 

76% 

83% 

— 


TABLE I. Comparing the hardest jet found in Z plus jets production at the LHC for R = 0.4. Shown is the percentage of jets 
found by the column algorithm that match the row algorithm at the 3% level, as measured by (1341) . 



J Ej 

XCONE N = 1 

SISCone-PR 

anti-fct 

C/A 

kt 

J E'J 

— 

64% 

72% 

64% 

45% 

42% 

XCONE N = 1 

63% 

— 

84% 

91% 

57% 

48% 

SISCone-PR 

71% 

85% 

— 

87% 

53% 

45% 

anti-fct 

63% 

91% 

87% 

— 

56% 

47% 

C/A 

44% 

56% 

53% 

56% 

— 

67% 

kt 

41% 

48% 

45% 

47% 

67% 

— 


TABLE If. Same as Table HJ but for R = 0.8 where larger differences are expected. 


Alternatively, one can define a stable cone to have its jet 
axis aligned with the true jet three-momentum, which is 
the default behavior in SISCone-PR [Hill- For small 
radius jets, the distinction is small, but we know of no 
1-jettiness measure (or jet function) that yields exactly 
stable cones by this second definition. 

Despite the above discussion, the mismatch in (l32l) is a 
small effect at small i?, with corrections that scale like 
or rn?/p^. Moreover, these finite R corrections are often 
subdominant to other effects present in practical jet al¬ 
gorithm implementations, such as the difference between 
local or global optimization or the treatment of jets at 
forward rapidities. 

In typical situations, maximizing Jhad, minimizing 
and finding stable cones (by either definition) yield 
rather similar jets, which are also similar to anti-Zct jets 
M- This is shown in Table U for a Z plus jets sample 
from Pythia 8.209 [H-H^ at the 14 TeV LHC, looking 
at the hardest jet with px > 50 GeV and R — 0.4 for all 
algorithms. The jet function approach is implemented 
with the Jet algorithm Q, which finds jets ordered by 
(1501) but we take the hardest jet by px- The 1-jettiness 
approach is implemented with XCONE (l^.[T7j|. using m 
with /3 = 20 The stable cone approach is implemented 
with SISCone-PR [H [3, taking the hardest jet by 
Px- The anti-fct approach is implemented with Fast Jet 
[28l| . For comparison, we also show results for the Cam¬ 
bridge/Aachen (C/A) [29M3Tj| and kt algorithms. 

As a similarity metric between the row and column jet al¬ 


Note that (I31II differs from the default XCONE measure based on 
Lorentz dot products. We use iterative one-pass minimization, 
which finds a local T^^ad from IRC safe seeds. 


gorithms in Table U we look at the fraction of the scalar 
Px not contained in the shared constituents of the two 
algorithms: 

^^i^rowDcolumn (34) 

^2Crow 4*T2 

Around 95% of the time, the four cone-like algorithms 
yield similar jet kinematics at the 3% level, while the 
differences with C/A and kt are noticeably larger. Be¬ 
cause all of the algorithms considered are longitudinally 
boost invariant, the similarities/differences between algo¬ 
rithms are essentially independent of the jet rapidity. In 
Table n we perform the same comparison for R = 0.8, 
where the differences are much larger as expected. 

Note that the study in Table |T] involves events that 
typically only have one hard jet without substructure. 
The four algorithms do exhibit larger variations when 
going to multiple jets, due to their differing treatment 
of jet overlap regions, though the algorithms still yield 
similar results when the jets are widely separated. For 
jets with substructure, interesting differences can arise 
for subjets at wide angles (see further discussion in i). 
Analogous to the discussion in Sec. lIVl two particles with 
equal momentum sharing will be clustered together if: 


^Ri2 iS 


R 

XCONE, anti-fct. 


y/2R 

J E'J' J 

(35) 

2R 

SISCone-PR, 



where 2 —> Rsep in the last line if there is sufficient smear¬ 
ing of the sub jet directions 0 As in Sec. lYl we can make 


^ The fact that the Je^ behavior is intermediate between anti-kt 
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SISCone-PR behave more like the Jet algorithm by 
switching from pT’-ordered to Jhad-ordered cones, and we 
can make XCONE behave more like the Jet algorithm by 
performing brute force global 7i minimization. In cases 
like this where the four algorithms give different results, 
the desired behavior depends on the physics of interest. 

VI. CONCLUSIONS 

The idea of jet finding as an optimization problem has 
a long history (see also [s^ - ldTj] !. so we find it remark¬ 
able that three seemingly different optimization strate¬ 
gies with different underlying philosophies correspond to 
the same meta optimization problem. As we saw, the ex¬ 
act correspondence only holds for fixed-angle cone jets, 
whereas at the LHC, complications arise when projecting 
four-vectors to the rapidity/azimuth plane. That said, 
the differences between the algorithms are small for small 
radius jets, and this correspondence helps explain why jet 
function maximization and 1-jettiness minimization yield 
such similar jets to stable cone jet algorithms. 

There are different conclusions one might draw from 
this analysis. To the extent that jet function maximiza¬ 
tion, 1-jettiness minimization, and stable cone finding 
yield similar jets (which are similar to anti-fct), one should 
simply choose the jet algorithm with the best computa¬ 


tional performance. Alternatively, to the extent that the 
three algorithms have different underlying philosophies, 
one might try to exploit those differences to develop new 
jet finding techniques with qualitatively different behav¬ 
iors or with better convergence in perturbative QCD cal¬ 
culations. Ultimately, one should choose the jet finding 
approach that best suits the intended physics application. 
In that spirit, this paper suggests that jet algorithms 
should be judged not by the elegance of their underlying 
philosophies but by the utility of their jet objects. 
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